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ABSTRACT 

The successful discovery of X-ray, optical and radio afterglows of gamma-ray bursts 
has significantly helped our understanding of these sources, and made possible the 
identification of host galaxies at cosmological distances. The energy release inferred 
in these outbursts rivals that of supernovae, while its photon energy output may 
considerably exceed it. Current models envisage this to be the outcome of a cataclysmic 
stellar event leading to a relativistically expanding fireball, in which particles are 
accelerated at shocks and produce nonthermal radiation. The substantial agreement 
between observations and the theoretical predictions of the fireball shock model 
provide confirmation of the basic aspects of this scenario. The continued observations 
show a diversity of behavior, providing valuable constraints for more detailed models. 
Crucial questions being now addressed are the beaming at different energies and its 
implications for the energetics, the time structure of the afterglow, its dependence on 
the central engine or progenitor system behavior, and the role of the environment on 
the evolution of the afterglow. 



1. Introduction 

The first discovery of a GRB afterglow at X-ray wavelengths was made with the BeppoSAX 



satellite (Costa et al, 1997) in February 1997, followed by a long-term detection and follow-up 



at optical wavelengths ( van Paradijs, et al., 1997| ). This represents a major and long-awaited 



breakthrough in the investigation of these mysterious sources, and emphasizes the importance 
of high resolution X-ray imaging techniques, which were succesfully used for other significant 
astronomical discoveries by Joachim Trumper and his colleagues with Rosat. The subsequent 
detection of other GRB afterglows followed in rapid succession, extending in some cases to radio 
and microwave wavelengths as well. This has made it possible to follow some of these sources 
over time scales of many months, making the identification of counterparts and host galaxies 
possible. The study of afterglows has provided strong confirmation for the generic fireball shock 
model of GRB, in which the 7-ray emission arises at radii of 10 13 — 10 16 cm ( [Rccs &: Meszaros 



1992] [Meszaros fc Rees, 1993| , [Rees k Meszaros , 1994) , |Paczyhski fe Xu, 1994 |Katz, 1994 |Sar 
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& Piran, 1995). In particular, this model led to the prediction of the quantitative nature of the 
signatures of afterglows ( Meszaros fc Rees, 1997a| ), in advance of the observational discoveries 
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1997) in GRB 970508 provided confirmation of the hypothesis that bursts were at cosmological 



distances. More recently, significant interest was aroused by the report ( |Kulkarni, et al, 199Sj ) of a 
spectroscopic redshift measurement of z = 3.4 for the afterglow of the burst GRB 971214, whose 
fluence would then correspond to a 7-ray energy of 10 53 ' 5 (f2 7 /47r) ergs, where AJ7 7 is the solid angle 
into which the gamma-rays are beamed. Such energies were discussed ( [Meszaros fc Rees, 1997b" ) in 
the context of compact mergers, such as neutron star-neutron star (NS-NS ) or black hole-neutron 
star (BH-NS ) mergers, which can power a relativistic fireball resulting in the observed radiation. 
In some of the detected afterglows there is evidence for a relatively dense gaseous environments, 
as suggested, e.g. by evidence for dust ( Reichart, 1998| ) in GRB970508, the absence of an optical 
afterglow and presence of strong soft X-ray absorption ( |Groot et al, 199~7| , Murakami, et al, 1997) 
in GRB 970828, the lack an an optical afterglow in the (radio-detected) afterglow (Taylor, et al.\ 



1997) of GRB980329, and spectral fits to the low energy portion of the X-ray afterglow of several 



bursts ( pwen, et al., 1998j ). The latter observations may be suggestive of "hypernova" models 
( Paczyhski , 1998 , Fryer fc Wooslcy, 1998| ), involving the collapse of a massive star or its merger 
with a compact companion, although it is probably too early to draw definite conclusions. While it 
is at present unclear which, if any, of these progenitors is responsible for GRB, or whether perhaps 
different progenitors represent different subclasses of GRB, there is general agreement that they 
all would be expected to lead to the generic fireball shock scenario mentioned above. Much of the 
current effort centers around trying to identify such progenitors more specifically, and trying to 
determine what effect, if any, they have on the observable fireball and afterglow characteristics. 



2. The Fireball Shock Scenario 



The isotropy of the angular distribution of GRB suggested already early on that GRB were 
outside our own galaxy ( Fishman &: Meegan, 1998| ), and the lack of structure associated with 
nearby galaxies indicated that their distance would be such that non-Euclidean and evolutionary 
effects may be important. This was also indicated by studies of the counts of the number of 
bursts per unit fluence or peak flux ( |Fishman &: Meegan, 19*98 , Horack, et al., 1996| ), as did 
evidence for time dilation in the bursts ([Norris, et al., 1995| ). Fits to cosmological counts including 
curvature, luminosity function and evolution effects indicated that the latter could be important 
( Reichart fc Meszaros , 1997| ). In particular, the evolution may be related to the star formation 
rate as a function of cosmological age ( [Totani, 1997 , Wijers et al., 199S , Krumholz, et al., 1998; ). 
The cosmological distances of some bursts have, since then (Mctzgcr et al, 1997, Kulkarni, etj 



al, 1998 , Djorgovski, et al., 1998 ) been directly measured. This, of course, poses significant 
constraints on the type of sources that may be able to produce this phenomenon. The typical 
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numbers characterizing the luminosity, duration, total energy per burst and event rate are 

L ~ 10 52 (S7/47r) erg s _1 , (1) 



n 



t w ~ 10 s , 
E~ 10 53 (fi/4vr) erg , 
KT 6 (n/47r) -1 galaxy'V 



(2) 
(3) 
(4) 



where VL is the solid angle into which the energy is channeled, and the durations t w have a large 
spread 10~ 3 ^ t w <, 10 3 s. The luminosities and total energies are also likely to be affected by a 
significant spread, whose magnitude is harder to estimate. The beaming solid angle Qj is also, 
thus far, very poorly constrained. 

An early suggestion for a GRB energy source at cosmological distances was the merger of 
NS-NS binaries (e.g. Paczyhski , 1986 , Goodman, 1986| ), whose estimated numbers and merger 
rates are close to the observed GRB detection rate ( |Narayan, Piran Shemi, 1991 , Phinney] 



1991 ). A related possibility are BH-NS binary mergers ( Narayan, Piran fe Shemi, 1991 , Paczyhski 



1992, Meszaros & Rees, 1992b), which have only slightly different event rates. An alternative 
suggestion ( IWoosley, 19931) was a "failed Supernova lb" , resulting from the collapse of fast rotating 
massive star which fails to produce a core collapse SN. Both a compact merger, and a fast rotating 
stellar collapse were known, from numerical simulations, to lead to a fast rotating torus around 
a central, high density object which eventually would develop into a black hole. The binding 
energy liberated in both of these events is of order 10 54 ergs, a good fraction of which will be 
carried away in a short pulse of vv and gravitational waves. If a small fraction of this emerges 
as electromagnetic energy, some of which is deposited in a region with sufficiently low baryon 
density, a relativistic fireball would form, whose radiation spectrum in the observer frame will 



peak in the MeV range ( Paczyhski , 1986 , Goodman, 1986 , Shemi fc Piran, 1990 ). One possible 
channel for converting some of this energy into electromagnetic form is given by the vv — > e + e~ 
process flEichler, et ai, 1989 ). A low baryon load condition (required to make the fireball highly 
relativistic) can occur naturally in a binary merger since a centrifugal barrier develops along the 
orbital symmetry axis, which is relatively free of baryons and provides an escape route for the 
e ± ,7 fireball flMcszaros fc Rees, 1992 ). This would also imply a collimation of the relativistic 
fireball, which would enhance the apparent flux (relative to what it would be if it were isotropic) 
by factors which conservatively might be of order 10-100. The same processes were estimated to 
be able to produce a collimated relativistic fireball in the failed SN lb model ( [Woosley, 1993j ). 
Another suggestion for powering a fireball was that this might result from magnetic flaring activity 
( Narayan, Paczyhski fc Piran, 1992; ) in the disrupted torus produced around the merged binary. 



Irrespective of the details of the progenitor, the resulting fireball is expected to be initially 
highly optically thick. From causality considerations the initial dimensions must be of order 
ct var ~> 10 7 cm, where t var is the variability timescale, and the luminosities must be much higher 
than a solar Eddington limit. Since most of the spectral energy is observed above 0.5 MeV, the 
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optical depth against 77 — > e ± is large, and an e , 7 fireball is expected. Due to the highly 
super-Eddington luminosity, this fireball must expand. Since in many bursts one observes a large 
fraction of the total energy at photon energies e 7 ^> lGeV, somehow the flow must be able to 
avoid degrading these photons ( since 77 — > would lead, in a stationary or slowly expanding 
flow, to photons just below 0.511 MeV [ Harding &: Baring, 1994| ). In order to avoid this, it seems 



inescapable that the flow must be expanding with a very high Lorentz factor, since in this case 
the relative angle at which the photons collide is less than T _1 and the threshold for the pair 
production is effectively diminished. Thus, photons with energy 

e 7 ,McV £5 104e t~MeV r 2 ■ ( 5 ) 

are able to escape, where T2 is the bulk Lorentz factor in units of 10 2 , and et is the energy 
of the target photons flMeszaros , 1995 ). Thus, simply from observations and general physical 



considerations, a relativistically expanding fireball is expected. However, the observed 7-ray 
spectrum is generally a broken power law, i.e., highly nonthermal. The optically thick e ± 7 fireball 
cannot, by itself, produce such a spectrum (it would tend rather to produce a modified blackbody). 
In addition, the expansion would lead to a conversion of internal energy into kinetic energy of 
expansion, so even after the fireball becomes optically thin, it would be highly inefficient, most of 
the energy being in the kinetic energy of the associated protons, rather than in photons. 

The most likely way to achieve a nonthermal spectrum in an energetically efficient manner is 
if the kinetic energy of the flow is re-converted into random energy via shocks, after the flow has 
become optically thin. This is a plausible scenario, in which two cases can be distinguished. In 
the first case (a) the expanding fireball runs into an external medium (the ISM, or a pre-ejected 
stellar wind flRees fc Meszaros , 1992| , Meszaros fc Rees, 1993| , [Katz, 1994| , [Sari & Piran, 1995| ). 



The second possibility (b) is that (even before such external shocks occur) internal shocks develop 



in the relativistic wind itself, faster portions of the flow catching up with the slower portions (Rees 



fc Meszaros , 1994 , Paczyhski fc Xu, 1994 ). In both of its versions, the fireball shock scenario is 
completely generic, independent of the specific nature of the progenitor, as long as it delivers the 
appropriate amount of energy (^ 10 52 erg) in a small enough region (^ 10 cm). This model has 
been successful in explaining the major observational properties of the gamma-ray emission, and 
is the main paradigm used for interpreting the GRB observations. 

External shocks will occur in an impulsive outflow of total energy E a in an external medium 
of average particle density n Q at a radius 



rde 



HFn^i/ 3^-2/3 cm ; (6) 



where the lab-frame energy of the swept-up external matter (T 2 m p c 2 per proton) equals the 
initial energy E a of the fireball, and rj = T = 10 2 ?72 is the final bulk Lorentz factor of the ejecta. 
The typical observer-frame dynamic time at the shock (assuming the cooling time is shorter 
than this) is tdec ~ fdec/cT 2 ~ seconds, for typical parameters, and the observed burst duration 
would be ~ tdec- The impulsive assumption requires that the initial energy input occur in a time 
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shorter than td ec - Variability on timescales shorter than tdec may occur on the cooling timescale 
or on the dynamic timescale for inhomogeneities in the external medium, but generally external 
shocks are not ideal for reproducing highly variable profiles ( pari Piran, 1998] ). However, they 
can reproduce bursts with several peaks ( [Panaitescu fc Meszaros , 1998a ) and may therefore be 



applicable to the class of long, smooth bursts, or bursts with a few smooth peaks. 

In a wind scenario, the energy input occurs non- impulsively over a timescale t w , and this can 
lead to internal shocks at smaller radii where i^ec < t w . One expects a similar initial dynamic 
behavior as in the impulsive injection case, assuming that a lab-frame luminosity L D and mass 
outflow M Q are injected at r ~ 77 and are continuously maintained over a time t w . Initially T oc r 
(with comoving temperature T oc r^ 1 ), followed by saturation at the value T max ~ rj which occurs 
at the radius r sat ~ /r/77, where r\ = L /M c 2 . In such wind models, internal shocks will occur at 
a radius ( [Rees Sz Meszaros , 1994| ) 

rdis ~ ct var r] 2 ~3x W u t var r]j cm, (7) 

where shells of different energies A77 ~ 7/ initially separated by ct var catch up with each other (with 
tvar < tw)- in order for internal shock radius r^is to occur above the wind photosphere radius 

r ph ~ M ( T T /(47rm pC r 2 ) ~ 10 14 L 53 % 3 cm, (8) 

but also at radii greater than the saturation radius r sa t ~ 777/ (so that most of the energy 
does not come out in a photospheric quasi-thermal radiation component) one needs to have 
7.5 x 10 1 L^tvaJ 5 ^ rj <, 3 x 10 2 L^tva/ 4 , where L53 = E^/t w . This type of models have the 



advantage (Rees &; Meszaros , 1994) that they allow an arbitrarily complicated light curve, the 
shortest variation timescale t var ^ 10 -3 s being limited only by the dynamic timescale at 77, where 
the energy input may be expected to vary chaotically. Such internal shocks have been shown 
explicitly to reproduce (and indeed even be required by) some of the more complicated light curves 



(Sari fc Piran, 1998], Panaitescu &: Meszaros , 1998d|; see however Dermer & Mitman, 1998). 



3. Progenitors and Central Engines 



Even before the measurement of a high redshift in a GRB afterglow, the difficulty in detecting 
the host galaxies of bright bursts(e.g, Schaefer, 1997 ) motivated the exploration of ways of 
increasing the possible total energy budget of GRB. The first explicit model to do this ( [Meszaros" 



Rees, 1997b| ), which was formulated before afterglows and their redshifts were discovered, involved 
converting a large fraction of the binding energy of a black hole and torus system (~ 10 54 ergs) 
into a fireball through MHD torques which power a Poynting jet outflow. Such a system would 
naturally arise from a NS-NS or a BH-NS merger (and it may also arise naturally in a failed 
SN lb or a "hypernova" ) . In the last year, a number of other possible energy sources have been 



considered as possible candidates for powering GRB (Fryer & Woosley, 1998, Popham, Woosley Sz 



Fryer, 1998). A fact which is not widely realized is that all plausible GRB progenitors suggested 
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so far (e.g. NS-NS or NS-BH mergers, Helium core - black hole [He/BH] or white dwarf - black 
hole [WD-BH] mergers, and the wide category labeled as hypernova or collapsars including failed 
supernova lb [SNe lb], single or binary Wolf-Rayet [WR] collapse, etc.) are expected to lead to a 
BH plus debris torus system, and they are all capable of producing relativistic outflows through 
the same mechanisms. An important point is that the overall energetics from these various 
progenitors do not differ by more than about one order of magnitude ( [Meszaros , Rees fc Wije~rs^ 



1998b|) 



Two large reservoirs of energy are available in a generic merger or collapse scenario: the 



binding energy of the orbiting debris, and the spin energy of the black hole flMeszaros Sz Rees 



1997b] , [Paczyriski , 1998]) 



The first mechanism, relying on the binding energy of the torus, can provide up to 42% of 
the rest mass energy of the disk, for a maximally rotating black hole. The vv — > e + e~ process 
can tap the thermal energy of the torus produced by viscous dissipation. For this mechanism to 
be efficient, the neutrinos must escape before being advected into the hole; on the other hand, the 
efficiency of conversion into pairs (which scales with the square of the neutrino density) is low if the 
neutrino production is too gradual. Typical estimates suggest a neutrino powered fireball of £ 10 51 
erg ( [Ruffert, et ai, 1997 , Popham, Woosley &: Fryer, 1998), except perhaps in the "collapsar" or 
failed SN lb case where estimates ( Popham, Woosley &: Fryer, 1998 ) indicate up to 10 52 ' 3 ergs for 
optimum parameters. If the fireball is collimated into a solid angle Vtj then of course the apparent 
"isotropized" energy would be larger by a factor (47r / ) , but unless Uj is ^ 10~ 2 — 10~ 3 this 
may fail to satisfy the apparent isotropized energy of 10 53 5 ergs implied by a redshift z = 3.4 
for GRB 971214 (a similar energy budget is implied by a possible photometric redshift z ~ 5 for 



GRB 980329, Fruchter, 1998j ). An alternative way to tap the torus energy is through dissipation 
of magnetic fields generated by the differential rotation in the torus ( [Paczyriski , 1992 , Narayan, 



Paczyriski & Piran, 1992, Meszaros & Rees, 1997b, Katz, 1997). Even before the BH forms, a 



NS-NS merging system might lead to winding up of the fields and dissipation in the last stages 
before the merger (Meszaros fc Rees, 1992], |Vietri, 1997a]). 



The second mechanism relies on tapping the spin energy of the black hole itself, and it can 
in principle extract up to 29% of the rest mass energy of the black hole. A hole formed from a 
coalescing compact binary is guaranteed to be rapidly spinning, and, being more massive, could 
contain more energy than the torus; the energy extractable in principle through MHD coupling to 



the rotation of the hole by the [Blandford Sz Znajek, 1977| effect could then be even larger than 
that contained in the orbiting debris ( Meszaros & Rees, 1997b] , paczyhski , 1998; ). Collectively, 
any such MHD outflows have been referred to as Poynting jets. 

The various progenitors (NS-NS , BH-NS , failed SNe, hypernovae and various collapsars) 
differ only slightly in the mass of the BH and that of the debris torus they produce, and they may 
differ more markedly in the amount of rotational energy contained in the BH. Strong magnetic 
fields, of order 10 15 G, are needed needed to carry away the rotational or gravitational energy in a 
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time scale of tens of seconds ( |Usov, l~994l , [Thompson, T994|) . 



If the magnetic fields do not thread the BH, then a Poynting outflow can at most carry the 
gravitational binding energy of the torus. For a maximally rotating and for a non-rotating BH this 
is 0.42 and 0.06 of the torus rest mass, respectively. The torus or disk mass in a NS-NS merger is 
( [Ruffert, et al, 1997Q M d ~ 0.1M© , and for a NS-BH, a He-BH, WD-BH merger or a binary WR 



collapse it may be estimated at ( Paczyriski , 1998| , [Fryer fe Woosley, 199S| ) M d ~ 1M . In the 



HeWD-BH merger and WR collapse the mass of the disk is uncertain due to lack of calculations 
on continued accretion from the envelope, so 1M@ is just a rough estimate. The energy available 
is then 

^ f 8 x lO 53 e(M d /M ) ergs , (fast rot.); 
J ' '' ' 1.2 X lO 53 e(M d /M ) ergs , (slow roi 

where e is the efficiency for converting gravitational into MHD jet energy. The largest energy 
reservoir is therefore likely to be associated with NS-BH, HeWD-BH or binary WR collapse, 
which have larger disks and fast rotation, the maximum energy being ~ 8 x lO 53 e(Md/M ) ergs; 
for the (fast rotating but smaller disk) NS-NS merger it is ~ 8 x lO 52 e(Md/O.lM ) ergs; and for 
the failed SNe lb (which is a slow rotator) it is ~ 1.2 x lO 53 e(Mrf/M ) ergs, Conditions for the 
efficient escape of a high-T jet may, however, be less propitious if the "engine" is surrounded by 
an extensive envelope, which is the case in the failed SNE lb or hypernova models. 

If the magnetic fields in the torus thread the BH, the spin energy of the BH can be extracted 
via the Blandford Znajek, 1977| (B-Z) mechanism (Meszaros & Rees, 1997b). The extractable 



energy is 

E ~ ef(a)M bh c 2 , (10) 

where e is the MHD efficiency factor and a = Jc/GM 2 is the rotation parameter, which equals 1 
for a maximally rotating black hole. The rotational parameter f(a) = 1 — ([1 + Vl — a 2 ]/2) 1 / 2 
is small unless a is close to 1, where it sharply rises to its maximum value /(l) = 0.29, so the 
main requirement is a rapidly rotating black hole, a <; 0.5. For a maximally rotating BH, the 
extractable energy is therefore 

Emax,bh ~ 0.29eM bh c 2 ~5x 10 53 e{M bh /M & ) ergs. (11) 

Rapid rotation is essentially guaranteed in a NS-NS merger, since the radius (especially for a soft 
equation of state) is close to that of a black hole and the final orbital spin period is close to the 
required maximal spin rotation period. Since the central BH will have a mass ( [Ruffert, et al., 



1997, Ruffert &: Janka, 1998| ) of about 2.5M , the NS-NS system can thus power a jet of up to 

Ema X ,bh ~ 1.3 x 10 54 e(M fch /2.5M Q ) ergs . (12) 

A maximal rotation rate may also be possible in a He-BH merger, depending on what fraction of 
the He core gets accreted along the rotation axis as opposed to along the equator ( [Fryer fc Woosley", 



1998), and the same should apply to the binary fast-rotating WR scenario, which probably does 
not differ much in its final details from the He-BH merger. For a fast rotating BH of 2.5 — 3M 
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threaded by the magnetic field, the maximal energy carried out by the jet is then similar or 
somewhat larger than in equation (12). The scenarios less likely to produce a fast rotating BH are 



the NS-BH merger (where the rotation parameter could be limited to a < Mng/M^, unless the 
BH is already fast-rotating) and the failed SNe lb (where the last material to fall in would have 
maximum angular momentum, but the material that was initially close to the hole has less angular 
momentum). The electromagnetic energy extraction from the BH in these could be limited by the 
f(a) factor, but a lower limit would be given by the energy available from the gravitational energy 
of the disk, in the second line of equation (^). 

Thus in the case of a Poynting (MHD) jet powered by the binding energy of the torus, the 
total energetics between the various models differs at most by a factor 20, whereas for Poynting 
jets powered by the spin energy of the black hole they differ by at most a factor of a few, depending 
on the rotation parameter. For instance, even allowing for low total efficiency (say 30%), a NS-NS 
merger whose jet is powered by the torus binding energy would only require a modest beaming of 
the 7-rays by a factor (An/Qj) ~ 20, or no beaming if the jet is powered by the B-Z mechanism, 
to produce the equivalent of an isotropic energy of 10 53 5 ergs. The beaming requirements of 
BH-NS and some of the other progenitor scenarios are even less constraining. Thus, even extreme 



redshifts z ~ 5 such as inferred by [Fruchter, 1998; can be easily satisfied by secanrios leading to a 
BH plus torus system. 

An interesting case is the apparent coincidence of GRB 980425 with the SN Ib/Ic 1998bw 
( Galama, et al, 1998| ). A simple but radical interpretation ( |Wang Wheeler, 1998]) is that 



all GRB may be associated with SNe Ib/Ic and differences arise only from different viewing 
angles relative to a very narrow jet. The difficulties with this are that it would require extreme 
collimations by factors 10 -3 — 10~ 4 , and that the statistical association of any subgroup of GRB 



with SNe Ib/Ic (or any other class of objects, for that matter) is so far not significant (Kippen, 



et al, 1998| ). If however the GRB 980425/ 1998bw association is real flWoosley, Eastman fc 



Schmidt, 1998| ), then we may be in the presence of a new subclass of GRB with lower energy 



E~f ~ 10 48 (J7j/47r) erg, which is only rarely observable even though its comoving volume density 
could be substantial. In this, more likely interpretation, the great majority of the observed 
GRB would have the energies E 1 ~ 10 54 (%/4tt) ergs as inferred from high redshift observations. 
Nonetheless, until further examples of such associations are discovered, one may not entirely 
discount the possibility of this being an extremely rare, low probability chance coincidence. 



4. Afterglows: the Simple "Standard" Model 

One can understand the dynamics of the afterglows of gamma-ray bursts in a fairly simple 
manner, without worrying about the uncertainties of the progenitor. This can be done through 
a relativistic generalization of the method used to understand supernova remnants without fully 
understanding the initiating explosion. The simplest hypothesis is that the afterglow is due to a 
relativistic expanding blast wave, which decelerates as time goes on ( [Meszaros fc Rees, 1997a ; 
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earlier simplified discussions were given by Katz, 1994b, Paczynski & Rhoads, 1993, Rees & 



Meszaros , 1992| ). The complex time structure of some bursts suggests that the central trigger 



may continue for up to 100 seconds. However, at much later times all memory of the initial time 
structure would be lost: essentially all that matters is how much energy and momentum has 
been injected; the injection can be regarded as instantaneous in the context of the much longer 
afterglow. Detailed calculations and predictions from such a model ( [Meszaros fc Rees, 1997a ) 



preceded the observations of the first afterglow detected, GRB970228 flCosta et al, 1997 



van 



Paradijs, et al, 1997j ) 



The simplest spherical afterglow model consist of a three-pieced power law spectrum with two 
breaks. At low frequencies there is a steeply rising synchrotron self-absorbed spectrum up to a 
self-absorption break u a , followed by a +1/3 energy index spectrum up to the synchrotron break 
v m corresponding to the minimum energy 7 m of the power-law accelerated electrons, and then a 
— (p — l)/2 energy spectrum above this break, for electrons in the adiabatic regime (where j~ p 
is the electron energy distribution above j m ). In addition, a third break is expected at energies 
where the electron cooling time becomes short compared to the expansion time, with a spectral 
slope — p/2 above that. With the inclusion of this third, "cooling" break 1/5, first calculated in 



Meszaros , Rees &: Wijers, 1998; and more explicitly detailed in pari, Piran fc Narayan, 1998| , one 
has what has come to be called the simple "standard" model of a GRB afterglow. This implicitly 
assumes spherical symmetry (for a jet with opening angle 9j this remains valid as long as T <; 0J ) 
and an impulsive energy input. As the remnant expands the spectrum simply moves to lower 
frequencies, so that in a given band the flux decays as a power law in time, whose index can 
change as breaks move through the observing band. 

The simple standard model has been remarkably successful at explaining the gross features 
of the GRB 970228, GRB 970508 and other afterglows ( [Wijers, Rees fc Meszaros , 1997| , [Tavani 



1997, Waxman, 1997] , Reichart, 1997). The multi-wavelength data analysis has in fact advanced 



to the point where one can use observed light curves at different times to extrapolate in time to 
get spectral snapshots at a fixed time (Waxman, 1997, Wijers fc Galama, 1998] ), allowing fits for 



the different physical parameters of the burst and environment, such as the total energy E, the 
magnetic and electron-proton coupling parameters eb and e e and the external density n Q . This has 
led to the temptation to take the assumed sphericity for granted. For instance, the lack of a break 
in the late light curve of GRB 970508 prompted the inference ( |Ramaprakash, A, et al, 1998 ) 



that all afterglows are essentially isotropic, leading to the very large (isotropic) energy estimate 
of 10 53 - 5 ergs in GRB 971214. However, what these fits constrain is only the energy per unit 
solid angle £ = (E/Qj). Furthermore, the assumption that a break should have been seen after 
months due to the fireball having become nonrelativistic, or T having dropped below an inverse jet 
opening angle size, is based upon the simple impusive energy input (a delta, or top hat function 
in initial energy and T, which are independent of angle). These are useful simplifications, but it is 
easy to see that departures from it are natural and would certainly not be surprising. Thus, as we 
emphasize below, there are so far no strong constraints on the possible anisotropy of the outflow 
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at various energies. 



5. Realistic Afterglows Models 

In the simplest departure from a spherical model the blast wave energy may be channeled 
into a solid angle Q,-. In this case one expects (Rhoads, 1997a, 1997b) a faster decay of T after 
sideways expansion sets in, and a decrease in the brightness is expected after the edges of the jet 

— 1/2 

become visible, when T drops below f2 ■ .A simple calculation using the usual scaling laws leads 
indeed to a steepening of the flux power law in time. The lack of such an observed downturn in 
the optical can, in a first approximation, be interpreted as an indication of the sphericity of the 
late stages of the fireball, and by extrapolation to the entire fireball including its early, gamma-ray 
emitting stages. 

There are, however, several important caveats. The first one is that the above argument 
assumes a simple, impulsive energy input (lasting ^ than the observed 7-ray pulse duration), 
characterized by a single energy and bulk Lorentz factor value. Estimates for the time needed 
to reach the non-relativistic regime, or T < Qj ^ few, could then be under a month ( [VietriJ 
1997a[) , especially if an initial radiative regime with T oc r -3 prevails. It is unclear whether, even 



when electron radiative time scales are shorter than the expansion time, such a regime applies, as 
it would require strong electron-proton coupling ([Meszaros , Rees &: Wijers, 1998| ). 

Furthermore, even the simplest reasonable departures from a top-hat approximation (e.g. 
having more energy emitted with lower Lorentz factors at later times, which still do not exceed 
the gamma-ray pulse duration) would drastically extend the afterglow lifetime in the relativistic 
regime, by providing a late "energy refreshment" to the blast wave on time scales comparable 

l_ j /■-) 

to the afterglow time scale ( [Rees Sz Meszaros , 1998[) . The transition to the r < O regime 



occurring at T ~ few could then occur as late as six months to more than a year after the outburst, 
depending on details of the brief energy input. Even in a simple top-hat model, more detailed 
calculations show that the transition to the non-relativistic regime is very gradual (5t/t <; 2) in 
the light curve. A numercial computation of the sideways expansion effects also shows that its 
effects are not so drastic as inferred from simple scaling for the material along the line of sight. 
This is because even though the flux from the head-on part of the remnant decreases faster, this is 
more than compensated by the increased emission measure from sweeping up external matter over 
a larger angle, and by the fact that the extra radiation, which arises at larger angles, arrives later 
and re-fills the steeper light curve. The sideways expansion thus actually can slow down the flux 



decay ( Panaitescu fc Meszaros , 1998c ) rather than making for a faster decay. Thus the conclusion 
is that we do not yet have significant evidence for whether the outflow is jet-like. The lack of a 
noticeable downturn in the light-curve is therefore, so far, compatible with either a spherical or a 
jet-like outflow. 



The ratio L^/L op t (or L^/L x ) can be quite different from burst to burst. The fit of Wijers 
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fe Galama, 1998 for GRB 970508 indicates an afterglow (X-ray energies or softer) energy per 
solid angle £52 = 3.7, while at z = 0.835 with /170 = 1 the corresponding 7-ray £527 = 0.63. On 
the other hand for GRB 971214, at z = 3.4, the numbers are £52 = 0.68 and £527 = 20. The 
7-ray bursts themselves require ejecta with T > 100. The gamma-rays we receive come only from 
material whose motion is directed within one degree of our line of sight. They therefore provide no 
information about the ejecta in other directions: the outflow could be isotropic, or concentrated in 
a cone of angle (say) 20 degrees (provided that the line of sight lay inside the cone). At observer 
times of more than a week, the blast wave would be decelerated to a moderate Lorentz factor, 
irrespective of the initial value. The beaming and aberration effects are less extreme so we observe 
afterglow emission not just from material moving almost directly towards us, but from a wider 
range of angles. 

The afterglow is thus a probe for the geometry of the ejecta — at late stages, if the outflow is 
beamed, we expect a spherically-symmetric assumption to be inadequate; the deviations from the 
predictions of such a model would then tell us about the ejection in directions away from our line 
of sight. It is quite possible, for instance, that there is relativistic outflow with lower T (heavier 
loading of baryons) in other directions; this slower matter could even carry most of the energy 
( Wijcrs, Rees fe Meszaros , 1997 , Paczyhski , 1998|) . This hypothesis is, in fact, supported to some 



degree by the fits of Wijers fe Galama, 1998 mentioned above 



6. Environment Effects, Evolution and Diversity 

One expects afterglows to show a significant amount of diversity. This is expected both 
because of a possible spread in the total energies (or energies per solid angle as seen by a given 
observer), a possible spread or changes in the injected bulk Lorenzt factors, and also from the fact 
that GRB may be going off in very different environments. 

The angular dependence of the outflow, and the radial dependence of the density of the 
external environment can have a marked effect on the time dependence of the observable afterglow 
quantities ( [Meszaros , Rees fe Wijers, 1998" ). So do any changes of the bulk Lorentz factor and 



energy output during even a brief energy release episode ( Rees fe Meszaros , 1998 ). The afterglow 
light curves are also affected by the degree of coupling between electrons and protons in the 
outflow (Meszaros , Rees & Wijers, 1998, Panaitescu fe Meszaros , 1998b| ). Diversity in the light 



curves of objects such as GRB 970508, e.g. sharp rises or humps followed by a renewed decay 
( Pcdersen, et al, 1998| , Piro, et al, 1998] ) is symptomatic of departures from the simple standard 



model. Detailed time-dependent model fits (panaitescu, Meszaros fe Rees, 1998 ) to the X-ray. 



optical and radio light curves of GRB 970228 and GRB 970508 show that, in order to explain the 
humps, a non-uniform injection or an anisotropic outflow is required. These fits indicate that the 
shock physics may be a function of the shock strength, and also indicate that dust absorption 
may be needed to simultaneously fit the X-ray and optical fluxes (the latter being affected more 
severely) . The effects of beaming (outflow within a limited range of solid angles) can be significant 
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(Panaitescu & Meszaros , 1998c), but are coupled with other effects, and a careful analysis is 
needed to disentangle them. 

Spectral signatures, such as atomic edges and lines, may be expected both from the outflowing 
ejecta (Meszaros & Rees, 1998a) and from the external medium (Perna & Loeb, 1998, Meszaros & 
Rees, 1998b ) in the X-ray and optical spectrum of afterglows. These may be used as diagnostics 
for the outflow Lorentz factor, or as alternative measures of the GRB redshift. An interesting 
prediction ( Meszaros fe Rees, 1998b| ; see also Ghisellini, et al., 1998 , Bottcher, et al., 1998|) is that 
the presence of a measurable Fe K-a emission line could be a diagnostic of a hypernova, since in 
this case one can expect a massive envelope at a radius comparable to a light-day where tt & 1, 
capable of reprocessing the X-ray continuum by recombination and fluorescence. 

The location of the afterglow relative to the host galaxy center can provide clues both for the 
nature of the progenitor and for the external density encountered by the fireball. A hypernova 
model would be expected to occur inside a galaxy, in fact inside a high density (n > 10 3 — 10 5 ). 
Some bursts are definitely inside the projected image of the host galaxy, and some also show 
evidence for a dense medium at least in front of the afterglow ( Pwen, et al., 19"98| ) . On the other 
hand, for a number of bursts there are strong constraints from the lack of a detectable, even faint, 
host galaxy ( |5chaefer, 1998 ). In NS-NS mergers one would expect a BH plus debris torus system 
and roughly the same total energy as in a hypernova model, but the mean distance traveled from 



birth is of order several Kpc ( Bloom, Sigurdsson fe Pols, 1998 ), leading to a burst presumably in 



a less dense environment. The fits of Wijers fe Galama, 1998 to the observational data on GRB 
970508 and GRB 971214 in fact suggest external densities in the range of n = 0.04-0.4 cm -1 , 



which would be more typical of a tenuous interstellar medium (however, Reichart 

- 3 ) 



Lamb, 1998 



report a fit for GRB 980329 with n Q ~ 10 cm A ). These could arise within the volume of the 
galaxy, but on average one would expect as many GRB inside as outside. This is based on an 



estimate mean NS-NS merger time of 10 years; other estimated merger times (e.g. 10 years, van 



den Heuvel, 1992] ) would give a burst much closer to the birth site. BH-NS mergers would also 
occur in timescales ^ 10 7 years, and would be expected to give bursts well inside the host galaxy 
dBloom, Sigurdsson fe Pols, 1998|) . 



7. Conclusions 

The simple blast wave model has proved quite robust in providing a consistent overall 
interpretation of the major features of gamma-ray bursts and their afterglows at various 
frequencies. However, the constraints on the angle-integrated energy, especially in 7-rays, are not 
strong, and beaming effects remain uncertain. Some of the observed light curve humps are likely 
to be indicative of either non-uniform injection episodes or anisotropic outflows. A relatively 
brief (1-100 s), probably modulated energy input appears the likeliest interpretation for most 
bursts, although in some progenitor scenarios there may be delayed effects. This can provide an 
explanation both for highly variable 7-ray light curves and late glitches in the afterglow decays. 
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One needs to be mindful of the possibility of there being more subclasses of classical GRB than 
just short ones and long ones. For instance, GRB with no high energy pulses (NHE) appear to 
have a different (but still isotropic) spatial distribution than those with high energy (HE) pulses 
(Pendleton, et al, 1996| ). Some caution is needed in interpreting this, since selection effects could 



lead to a bias against detecting HE emission in dim bursts (Norris, 1998). A connection to peculiar 



supernova events, or the formation of a temporarily rotationally stabilized strong-field pulsar are 
also possibilities which cannot be ruled out, 

Significant progress has been made in understanding how gamma-rays can arise in fireballs 
produced by brief events depositing a large amount of energy in a small volume, and in deriving 
the generic properties of the long wavelength afterglows that follow from this. There still remain 
a number of mysteries, especially concerning the identity of their progenitors, the nature of the 
triggering mechanism, the transport of the energy and the time scales involved. Nevertheless, 
even if we do not yet understand the intrinsic gamma-ray burst central engine, they may be the 
most powerful beacons for probing the high redshift (z > 5) universe. Even if their total energy 
is reduced by beaming to a "modest" ~ 10 52 — 10 52 ' 5 ergs in photons, they are the most extreme 
phenomena that we know about in high energy astrophysics. The modeling of the burst mechanism 
itself will continue to be a major challenge to theorists and to computational techniques. However, 
there is every prospect for continued and vigorous developments both in the observational analyses 
and in the theoretical understanding of these fascinating objects. 

I am grateful to Martin Rees for stimulating collaborations on this subject, as well as to 
Ralph Wijers, Hara Papathanassiou and Alin Panaitescu. This research is supported in part by 
NASA NAG5-2857 
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